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Abstract 

We study numerically the maximal and minimal height distributions (MAHD, MIHD) of the 
nonlinear interface growth equations of second and fourth order and of related lattice models 
in two dimensions. MAHD and MIHD are different due to the asymmetry of the local height 
distribution, so that, in each class, the sign of the relevant nonlinear term determines which one 
of two universal curves is the MAHD and the MIHD. The average maximal and minimal heights 
scale as the average roughness, in contrast to Edwards- Wilkinson (EW) growth. All extreme 
height distributions, including the EW ones, have tails that cannot be fit by generalized Gumbel 
distributions. 
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Extreme value statistics (EVS) has already been applied in several fields of science and 

n n n 

engineering pj, 12|, |3|]. It has recent important applications in surface science, e. g. for 
modeling the evolution of corrosion damage at time scales not easily accessible to experiment 
^. In uncorrelated random variable sets, the statistics of the nth extrema is described by 
the Gumbel's first asymptotic distribution if the probability density functions (PDF) 
of those sets decrease faster than a power law. However, deviations from this statistics are 
expected in fluctuating interfaces if there are strong correlation of local heights. In one 
dimension, this is the case of Edwards- Wilkinson (EW) interfaces (Brownian curves) 
and other Gaussian interface models 8|. On the other hand, fluctuations of other global 
quantities in various physical systems follow generalized Gumbel distributions ^, [lo|, 
i. e. the first asymptotic distribution with noninteger n values. This is explained by the 



connections between the EVS o 
drawn from exponential PDF [12 



correlated variables and sums of independent variables 
which shows that Gumbel statistics goes far beyond the 
description of uncorrelated variables sets. 

The interface models where maximal height distributions (MAHD) were previously cal- 
culated [6|, 17|, 18| are symmetric with respect to the average height, consequently MAHD and 
minimum height distributions (MIHD) are the same. However, there is a large number of 
real interfaces where the up-down symmetry is broken 13 L such as those described by the 



nonlinear growth models of Kardar-Parisi-Zhang (KPZ) |l^ and of Villain-Lai-Das Sarma 
(VLDS) for molecular beam epitaxy 15(], which raises the question whether MAHD and 
MIHD are the same in those systems. This is particularly important for two-dimensional 



interfaces due to the variety of real growth processes which show KPZ |16l . 1171 ] and VLDS 
scaling [l^. Recent works on persistence in VLDS growth also motivate such study, since 
different exponents for positive and negative height persistence were obtained 1^. More- 
over, the above scenario raises additional (and not less important) questions for the nonlinear 
models. The first one is connected to the possibility of fitting their extreme height distri- 
butions (EHD) by generalized Gumbel distributions, similarly to other correlated systems. 
The second one is the scaling of the average maximal height, since EW interfaces showed 
an unanticipated scaling as the square of the average roughness [nl], in contrast to several 
one-dimensional interfaces. This is essential to correlate surface roughness with the extreme 
events. 



The aim of this letter is to address those questions by performing a numerical study of 
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the MAHD and MIHD in the steady states of the KPZ and VLDS equations and of various 
lattice models belonging to those classes in 2 + 1 dimensions. We will show that, for each 
growth class, two universal distributions are obtained, which may be a MAHD or a MIHD 
of a given model depending on the sign of the coefficient of the relevant nonlinear term. 
Combination of data collapse and extrapolation of amplitude ratios (e. g. skewness and 
kurtosis) of those distributions are used to separate systems with coefficients of different 
signs. In order to illustrate the drastic effects that asymmetric PDF (i. e. distributions of 
local heights) may have on MAHD and MIHD, we will discuss their differences in a random 
deposition-erosion model on an inert fiat substrate. We will also show that average maximal 
and minimal heights is all those models scale as the average roughness, as usually expected. 



which shows that the EW scaling is an exception 



Finally, we will show that KPZ, 



VLDS and EW distributions cannot be fit by generalized Gumbel distributions. 

MAHD and MIHD were calculated for the KPZ equation ^ = i/gV^/i + A2(V/i)^ + t]{x, t) 
[{r] (x, t) r]{x', t')) = D5'^{x - x')S {t - t')] in dimension d = 2, with z/2 = 0.25, D = 5x 10~^ 
and A2 = V75 {g = \2'D/i'2^ = 24), in discretized boxes with spatial step Ax = 1, time 
increment At = 0.04 and linear sizes 8 < L < 64. A simple Euler integration method [2^ 
and a scheme for suppression of instabilities [21] were adopted. We also simulated three 
discrete KPZ models in sizes 32 < L < 256: the restricted solid-on-solid (RSOS) model 



22l |. the ballistic deposition (BP) 23l | and the etching model of Mello et al l2J]. From 



inspection of their growth rules 25|] one knows that A2 > for BD and the etching model 
and A2 < for the RSOS model. The VLDS equation § = -v^V^h + X^V^iVhf + r]{x,t) 
was integrated with z/4 = 1, A4 = 1, = 1/2, and At = 0.01, using the same methods, in 
sizes 8 < L < 32. We also simulated a generalized conserved RSOS model (CRSOS) |26| . 
whose original version was proposed in Ref. 27|) and which belongs to the VLDS class, in 
sizes 16 < L < 128. The EW equation (KPZ with A2 = 0) was integrated with v = 1.5, 
D = 1/2 and At = 0.01 in box sizes 8 < L < 64. For each model and each lattice size, 
distributions with at least 10^ different configurations were obtained to ensure high accuracy, 
which is particularly important at their tails. The extremes were calculated relatively to 
the average height of each configuration, the minima being absolute values of the differences 
from the average. 

In Fig. la we show the scaled MAHD and MIHD of the KPZ equation in box size L = 64. 
In these plots, P{m)dm is the probability that the extreme lies in the range [m,m + dm], 
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X = {m — (m)) /a and cr = (^{m?) — {m)'^^ . The high accuracy in Fig. la allows us to 
distinguish those curves (log-linear plots also show discrepancies in the right tails). Results 
for smaller box sizes show that the finite-size effects are negligible, confirming that MAHD 
and MIDH are actually different. 

Gumbel's first asymptotic distribution used to compare our data is g (x; n) = 
uexp(^—n e~^^^~^'^'' + b {x + s) ^, where b = ^ip' (n), s = [Inn — '?^(n}l /6 and u = 
n"'b/r{n), with T {x) the Gamma function and {x) = dlnT {x)/dx 9|, [ill]. In Fig. lb, 
we show the MAHD of the KPZ equation and the Gumbel curve with the same skewness 
0.79, which has n = 1.95. Although the fit near the peak is reasonable, there are significant 
differences in the tails. Data for different lattice sizes in Fig. lb show that discrepancies 
are not consequence of finite-size effects. Similar disagreement is observed when we try to 
fit the MIHD with a Gumbel curve of skewness 0.65 [n = 2.75). However, the right tails of 
MAHD and MIHD of the KPZ class tend to simple exponential decays for large m, similarly 
to the Gumbel curves. 

In Fig. 2a we show the scaled MAHD of the integrated KPZ equation and of BD, and 
the MIHD of the RSOS model. In Fig. 2b we show the MIHD of the KPZ equation and of 
the etching model, and the MAHD of the RSOS model. There is excellent data collapse in 
both plots, which confirm that MAHD (MIHD) of models with A2 > are equal to MIHD 
(MAHD) of models with A2 < 0. This illustrates the possibility of using the MAHD and 
MIHD to identify the sign of the coefficient of the nonlinear term in cases where it is not 
known a priori. At this point, EVS is superior to the scaling of the local height distributions 
(the PDF), which might also reveal the sign of the nonlinear terms if distortion by huge 
finite-size effects were not so frequent (28| . 

The visual agreement between those distributions and the small finite-size effects are 
quantitatively confirmed by estimates of their skewnesses and kurtosis in various system 
sizes. Figs. 3a and 3b show the skewness of the same models of Figs. 2a and 2b, respectively, 
as a function of 1/ L^/^ (BD data were not shown in Figs. 3a and 3b because they superimpose 
the etching model data). The small finite-size dependence of the data for the KPZ equation, 
BD and the etching models leads to S" ^ 0.79 for (A2 > 0)-MAHD and (A2 < 0)-MIHD, and 
S ^ 0.65 for (A2 > 0)-MIHD and (A2 < 0)-MAHD. Surprisingly, the largest finite-size effects 
are observed in the RSOS data, which uses to be the best discrete KPZ model for numerical 
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study of roughness scaling 



The MAHD and MIHD for the VLDS equation in box size L = 32 are plotted in Fig. 4a. 
Differences in the peaks are tiny, but discrepancies in the tails are clearly observed. Finite- 
size effects are also negligible in this case. The MAHD has skewness S ~ 0.63 and the MIHD 
has S ~ 0.55. In contrast to the KPZ models, those distributions have Gaussian-shaped 
right tails [exp(— m^)]. For this reason, fits with generalized Gumbel distributions (which 
tend to simple exponentials as m ^ oo) are not possible. This is confirmed for MAHD in 
Fig. 4a by comparison with Gumbel's curve with n = 2.90, which has skewness 0.63. 

For the CRSOS model, MAHD and MIHD show significant finite-size effects, similarly 
to the RSOS model. However, extrapolation of the skewness of distributions in finite-size 
lattices, shown in Fig. 4b, suggest that MAHD and MIHD of both models are asymptotically 
the same. This means that A4 > for the CRSOS model, similarly to its one-dimensional 
original version, which is exactly solvable (29| . 

The difference between MAHD and MIHD can be easily explained in a model of random 
deposition and erosion with an inert flat substrate in the erosion-dominated regime. Let 
q > 1/2 be the probability of single-particle erosion and 1 — g of deposition, and assume that 
erosion is possible only if h > 0. A steady state is attained with average height (relatively 
to the substrate) {h) = q/{2q-l) and PDF P (h) oc exp{-h/{h)) (for q close to 1/2). 
Now consider that one measures the extremes in a set of L (independent) columns. With 
that PDF, MAHD is given by Gumbel's first asymptotic distribution with n = 1. For 
large L, the minimum absolute height is typically at the substrate, thus fluctuations of the 
relative minima are dominated by fluctuations of the average height, which are Gaussian 
[oc exp (^—m'^/ay/L^ , a = ^''^^ so it is the MIHD. The difference between MAHD 

and MIHD is easily conflrmed by visual inspection of the plots of these functions. It is 
related to the highly asymmetric local height distribution of this model (skewness of PDF 
is SpDF = 2), in contrast to the slight asymmetry of KPZ {Spdf ~ 0.26 |28|) and VLDS 

{spDF ~ 0.20 y). 

Now we analyze the average values of extremes of KPZ and VLDS interfaces. They are 
assumed to scale as (m) ~ L"™ , while the average roughness scales with the roughness expo- 
nent a. We estimate am by extrapolation of effective exponents (L) 



ln[(m)(L)/(m)(L/2)] 
In 2 ■ 



as shown in Figs. 5a (KPZ) and 5b (VLDS). The estimates of am are consistent with the best 
known estimates of the roughness exponents a ~ 0.39 (KPZ) 28|] and a ~ 0.67 (VLDS) 
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261]. Consequently, the average values of extremes scale as the average roughness in the 



KPZ and VLDS classes in 2 + 1 dimensions. This contrasts with the scaling as the squared 
roughness in the EW class [11], which we also confirmed by simulation. 

Finally, in Fig. 6 we show the MAHD of the EW equation and the generalized Gumbel 
distribution with n = 2.6, which has the same skewness. Despite the good agreement in 
almost three decades of the scaled P{m), the discrepancy in the tails is clear. Again, data 
for two box sizes {L = 64 and L = 32) show that this is not caused by finite-size effects 
nor to low accuracy of the data. On the other hand, the analytical prediction by Lee 11 1 
of a Gaussian-shaped tail [~ exp (— m^)] of the MAHD is confirmed by the trend of our 
data for large m. Together with the above results for KPZ and VLDS classes, it shows that 
EVS of important interface growth models in two dimensions are not connected to the EVS 
of independent variables, despite the wide applicability of Gumbel statistics to correlated 
systems. 

In summary, we showed that interface growth models with asymmetric local height dis- 
tributions have different maximal and minimal height distributions, the most important 
examples being the KPZ and the VLDS classes in two dimensions. In each class, a pair of 
universal curves may be maximal of minimal height distributions depending on the sign of 
the relevant nonlinear term. The average maximal and minimal heights of KPZ and VLDS 
models scale as the average roughness, in contrast to the EW class. All extreme height 
distributions, including the EW ones, cannot be fit by generalized Gumbel distributions. 
Although most works on statistical properties of interfaces focus on features of height dis- 
tributions and/or roughness scaling [13|, recent studies show that the statistics of global 
quantities are very useful to characterize real growth processes {2, |3o|. The EVS has the 
same advantages of roughness distribution scaling for this task, such as weak finite-size ef- 
fects, and also reveals the sign of the nonlinear terms. Information on rare events is also 
essential in systems where drastic changes in the dynamics occur if the global minima or 
maxima attain certain values, such as in corrosion damage. On the other hand, for some 
applications (from friction to parallel computing) the distributions of local extremes may be 
important, and the present study certainly motivates additional studies of those quantities 
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32|. 



6 



Acknowledgments 



TJO acknowledges support from CNPq and FDAAR acknowledges support from CNPq 
and FAPERJ (Brazilian agencies). 



[1] E. J. Gumbel, Statistics of Extremes (Columbia University, New York, 1958). 
[2] S. T. Bramwell et al, Nature (London) 396, 552 (1998). 

[3] J. P. Bouchaud and M. Mezard, J. Phys. A 30, 7997 (1997); R. W. Katz et al. Adv. Water 
Resour. 25, 1287 (2002); S. T. Bramwell et al, Europhys. Lett. 57, 310 (2002); J. F. Eichner 
et al, Phys. Rev. E 73, 016130 (2006); S. Redner and M. R. Petersen, Phys. Rev. E 74, 061114 
(2006). 

[4] G. Engelhardt and D. D. Macdonald, Corros. Sci. 46, 2755 (2004). 

[5] R. A. Fisher and L. A. Tippett, Proc. Cambridge Philos. Soc. 24, 180 (1928). 

[6] S. Raychaudhuri et al, Phys. Rev. Lett. 87, 136101 (2001). 

[7] S. N. Majumdar and A. Comtet, Phys. Rev. Lett. 92, 225501 (2004); J. Stat. Phys. 119, 777 

(2005) . 

[8] G. Gyorgyi et al, Phys. Rev. E 68, 056116 (2003); G. Schehr and S. N. Majumdar, Phys. Rev. 
E 73, 056103 (2006); G. Gyorgyi et al, Phys. Rev. E 75, 021123 (2007). 

[9] S. T. Bramwell ct al, Phys. Rev. Lett. 84, 3744 (2000). 
[10] T. Antal et al, Phys. Rev. Lett. 87, 240601 (2001). 
[11] D.-S. Lee, Phys. Rev. Lett. 95, 150601 (2005). 

[12] E. Bertin, Phys. Rev. Lett. 95, 170601 (2005); E. Bertin and M. Clusel, J. Phys. A 39, 7607 

(2006) . 

[13] A.L. Barabasi and H.E. Stanley, Fractal concepts in surface growth (Cambridge University 

Press, Cambribge, England, 1995). 
[14] M. Kardar, G. Parisi and Y.-C. Zhang, Phys. Rev. Lett. 56, 889 (1986). 
[15] J. Villain, J. Phys. I 1, 19 (1991); Z.-W. Lai and S. Das Sarma, Phys. Rev. Lett. 66, 2348 

(1991). 

[16] J. Krim and G. Palasantzas, Int. J. Mod. Phys. B 9, 599 (1995). 

[17] R. Paniago et al, Phys. Rev. B 56, 13442 (1997); M. C. Salvadori et al, Phys. Rev. E 58, 6814 



7 



(1998); M. U. Kleinke et al, Appl. Phys. Lett. 74 1683 (1999); P. L. Schilardi et al, Phys. 

Rev. B 59, 4638 (1999); A. E. Lita and J. E. Sanchez, Jr., Phys. Rev. B 61 7692 (2000); D. 

Tsamouras et al, Appl. Phys. Lett. 79, 1801 (2001). 
[18] Y.-L. He et al, Phys. Rev. Lett. 69, 3770 (1992); C. Thompson et al, Phys. Rev. B 49, 4902 

(1994); D. C. Law et al, J. Appl. Phys. 88, 508 (2000). 
[19] M. Constantin et al, Phys. Rev. E 69, 061608 (2004). 
[20] K. Moser et al, Physica A 178, 215 (1991). 
[21] C. Dasgupta et al, Phys. Rev. E 55, 2235 (1997). 
[22] J. M. Kim and J. M. Kostcrlitz, Phys. Rev. Lett. 62, 2289 (1989). 
[23] M. J. Void, J. Coll. Sci. 14 (1959) 168; J. Phys. Chem. 63, 1608 (1959). 
[24] B. A. Mello et al, Phys. Rev. E 63, 041113 (2001). 
[25] W. E. Hagston and H. Ketterl, Phys. Rev. E 59, 2699 (1999). 

[26] F. D. A. Aarao Reis, Phys. Rev. E 70, 031607 (2004); F. D. A. Aarao Reis and D. F. Frances- 

chini, Phys. Rev. E 61, 3417 (2000). 
[27] Y. Kim et al, J. Phys. A: Math. Gen. 27, L533 (1994). 

[28] E. Marinari et al, J. Phys. A 33, 8181 (2000); F. D. A. Aarao Reis, Phys. Rev. E 69, 021610 
(2004). 

[29] S.-C. Park ct al, Phys. Rev. E 65, 015102(R) (2002). 
[30] S. Moulinet et al, Phys. Rev. E 69, 035103(R) (2004). 
[31] Z. Toroczkai et al, Phys. Rev. E 62, 276 (2000). 
[32] F. Hivert et al, J. Stat. Phys. 126, 243 (2007). 



8 



0.5 




-4 -2 2 4 6 8 

X 



FIG. 1: (Color online) a) Scaled MAHD (solid curve) and MIHD (dashed curve) of the KPZ 
equation in box size L = 64. b) Scaled MAHD of the KPZ equation (solid curve for L = 64, 
squares for L = 32) and generalized Gumbel distribution with n = 1.95 (dashed curve). 
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FIG. 2: (Color online) a) Scaled MAHD of the KPZ equation (solid curve) and BD (triangles), 
and MIHD of the RSOS model (squares), b) Scaled MIHD of the KPZ equation (solid curve) and 
etching model (triangles), and MAHD of the RSOS model (squares). For KPZ equation, box size 
is L = 64, and for discrete models L = 256. 
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FIG. 3: a) Finite-size dependence of the skewness S of MAHD of the KPZ equation (squares) and 
etching model (triangles), and of MIHD of the RSOS model (crosses), b) Finite-size dependence 
of S of MIHD of the KPZ equation (sqTiares) and etching model (triangles), and of MAHD of the 
RSOS model (crosses). In both plots, the variable in the abscissa was chosen to make clearer the 
evolution of the data as L — ^ oo. 



11 



-4 



-2 2 4 
X 



6 



0.4- ^ 
0.3- 




e 



-X- 
€)- 



b) 




0.2 



0.3 



1/L 



2/3 



FIG. 4: (Color online) a) Scaled MAHD (solid curve) and MIHD (dashed curve) of the VLDS 
equation in L = 32 and generalized Gumbel distribution with n = 2.90 (dotted curve), b) Finite- 
size dependence of the skewness S of EHD of the VLDS equation (crosses for MAHD, circles for 
MIHD) and of the CRSOS model (squares for MAHD, triangles for MIHD). The variable l/L^/^ 
provides the best linear fits of the data (dashed lines). 
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FIG. 5: Finite-size dependence of effective exponents Om {L) of discrete models: a) KPZ class (BD: 
crosses; etching: triangles; RSOS: squares) and b) VLDS class (CRSOS model). L is the average 
size among L and L/2. The variables in the abscissa provide the best linear fits (dashed lines) with 
exponents A = 1/2 (a) and A = 2/3 (b). 




FIG. 6: (Color online) Scaled MAHD of the EW equation in L = 64 (solid curve) and L = 32 
(squares) and the generalized Gumbel distribution with n = 2.6 (dashed curve). 
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